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Аннотация 

Введение. Математическое моделирование эффективно при анализе промышленной безопасности на 
металлургических предприятиях, в частности для отслеживания проблем системы «человек — машина». Чтобы 
ввести в рассмотрение фактор времени, задействуют рекуррентные соотношения (в дискретной модели) и 
дифференциальные (в непрерывной). Однако необходимо также решить проблему привязки параметров модели 
к реальным условиям производственной среды и к человеческому фактору. Цель данного 
исследования — создание метода определения параметров имитационных математических моделей динамики 
психофизиологических показателей оператора, влияющих на его работу. 

Материалы и методы. Психофизиологическое состояние (ПФС) оператора оценивали по работоспособности, 
утомляемости и ошибаемости. Данные собрали по тесту цифровой корректурной пробы (ЦКП). На основании 
полученных результатов вычислили экспериментальные значения показателей ПФС оператора, которые привели 
к нормированной шкале [0, 1]. Эти показатели для конкретного респондента, математическую модель и 
разработанный алгоритм задействовали при определении числовых значений параметров модели. Для 
интерпретации показателей работоспособности, утомляемости и ошибаемости ввели шкалы с пятью градациями. 
Результаты исследования. Использование модифицированного авторами варианта математической модели 
показало значительное улучшение ее прогностических свойств. Из 10 участников наилучший результат оказался у 
респондента № 7, худший — у респондента № 8. В течение 1-го часа работы (с 9.00 до 10.00) их работоспособность 
выросла примерно одинаково, с 0,5—0,55 почти до 0,6. Затем показатель респондента № 7 активно увеличивался и до 
конца рабочего дня оставался существенно выше уровня «хороший». Показатель респондента № 8 падал и с 14.00 до 
15.00 оказался ниже среднего. Разницу во многом определили хронотипы операторов. Их хронофизиологические 
особенности сказались также на утомляемости и ошибаемости. Для разных участников экспериментов 
варьировалось качество модели. В одном случае оно оказалось отличным (средняя относительная ошибка <5 %), в 
трех случаях — хорошим (<10 %), в четырех — удовлетворительным (<15 96). 

Обсуждение и заключение. Предлагаемый подход дает возможность получить для каждого индивидуума 
динамические профили его психофизиологических характеристик, оценить их взаимосвязи и выполнить прогноз 
на основе модифицированной математической модели. Однако для расширения функциональных возможностей 
моделей в реальных условиях работы оператора металлургического оборудования требуется увеличить объем 
выборки, уменьшить шаг дискретного времени и выполнить исследования для различных условий работы с 
учетом технологических, климатических, экологических, психологических и прочих факторов. 
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Abstract 

Introduction. Mathematical modelling is effective in the analysis of industrial safety at metallurgical plants, in 
particular for tracking problems of the man — machine system. To introduce the time factor, recurrence relations (in a 
discrete model) and differential relations (in a continuous model) are used. However, it is also necessary to solve the 
problem of linking the model parameters to the real conditions of the production environment and to the human factor. 
The aim of this study is to create a method for determining the parameters of simulation mathematical models of the 
dynamics of the operator's psychophysiological indicators affecting the work. 

Materials and Methods. The operator's psychophysiological state (PPS) was assessed by performance, fatigue levels, 
and error rate. The data were collected by the Digital Correction Task (DCT) test. Based on the obtained results, the 
experimental values of the operator's PPS indicators, which were reduced to the normalized scale [0, 1], were 
calculated. These indicators for a particular respondent, the mathematical model and the developed algorithm were used 
to determine the numerical values of the model parameters. In order to interpret the indicators of performance, fatigue 
and error rate, we introduced scales with five gradations. 

Results. The use of the authors" modified version of the mathematical model showed a significant improvement in its 
prognostic properties. Out of 10 participants the best result was shown by respondent no. 7, the worst result was shown 
by respondent no. 8. During the first working hour (from 9.00 to 10.00) their performance increased almost equally, 
from 0.5—0.55 to almost 0.6. Then the score of respondent no. 7 increased and remained well above the “good” level 
until the end of the day. The score of respondent no. 8 dropped and was below average from 14.00 to 15.00. The 
difference was largely determined by the operators’ chronotypes. Their chronophysiological characteristics also affected 
fatigue and error rate. The model's quality varied for different participants in the experiments. In one case it was 
excellent (mean relative error <5%), in three cases it was good (<10%) and in four it was satisfactory (<15%). 
Discussion and Conclusion. The proposed approach allows us to obtain the dynamic profiles of psychophysiological 
characteristics for every individual, to assess their interrelationships and to perform a prediction on the basis of a 
modified mathematical model. However, in order to extend the functionality of the models to the real working 
conditions of the metallurgical plant operator, it 1s necessary to increase the sample size, reduce the discrete time step 
and conduct studies for different working conditions, considering technological, climatic, environmental, psychological 
and other factors. 


Keywords: safety at metallurgical enterprises, man — machine system, psychophysiological state of the operator, 
operator chronotype, operator performance 
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Введение. Согласно данным Федеральной службы по экологическому, технологическому и атомному 
надзору Российской Федерации, к началу 2022 года на металлургических и коксохимических производствах 
зарегистрировано 1280 опасных производственных объектов, 19 из которых относятся к первому классу 
опасности, 325 — ко второму, и 922 — к третьему!. Как пишет журнал «Промышленные страницы»?, в 
2022 году на российских металлургических предприятиях самыми частыми происшествиями были пожары и 
взрывы. Эти инциденты обусловлены техническими причинами и нарушениями правил промышленной 
безопасности. Например, в январе 2022 года трое сотрудников Новолипецкого металлургического комбината 
проводили плановый осмотр межцехового газопровода без обязательной газозащитной аппаратуры. Рабочие 
погибли в результате отравления токсичными парами и газом. 

B [1] выделены технические, санитарно-гигиенические, организационные и психофизиологические группы 
факторов аварийности и травматизма на металлургических и коксохимических предприятиях. Последние две 
группы тесно связаны с человеком — участником производственного процесса. С точки зрения 
производственной безопасности металлургических предприятий особая ответственность возлагается на 
операторов высокотехнологичных агрегатов. В инженерной психологии их принято называть 
«человек-оператор». В рамках представленной работы в этом значении будем употреблять слово «оператор». 
Иногда именно он предотвращает переход опасности в инцидент или аварию на производстве. 

ГОСТ 12.0.003—2015? выделяет в самостоятельный блок опасные и вредные производственные факторы 
психофизиологического воздействия на человека. Оператор испытывает нервно-психические перегрузки, 
связанные с напряженностью трудового процесса. В их числе: 

— умственное перенапряжение; 

— перенапряжение анализаторов; 

— монотонность труда; 

— неуверенность в действиях из-за недостатка образования и опыта. 

Эмоциональные перегрузки ведут к переутомлению, плохому самочувствию, стрессовому состоянию и T. п. [1]. 

Надежность человека как элемента сложной технической системы зависит от внутренних и внешних 
условий, которые меняются во времени. Есть еще одна переменная — сам человек. Работу выполняют люди с 
разными личностными качествами, здоровьем, опытом и пр. 

В [2] исследовались сведения о несчастных случаях (НС) в кузнечно-прессовом производстве (КПП) за 
1985, 1987 и 1989 гг. Отмечены некоторые особенности в распределении по количеству часов от начала работы 
до инцидента. Эти закономерности можно объяснить человеческим фактором, точнее суточным ритмом 
трудовой деятельности (рис. 1). 
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Рис. 1. Распределение НС в кузнечно-прессовом производстве по количеству часов от начала работы до НС [2] 


! Годовой отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору в 2021 году. Москва: 2022. 
С. 407. URL: https://www.gosnadzor.ru/public/annual терогіѕ/Годовой%20отчет%20за%202021%20г.ра (дата обращения: 22.11.2023). 

2 Опасная тенденция: аварии и ЧС на металлургических производствах в 2022-2023 годах. Промышленные страницы. URL: 
https://indpages.ru/safe/avareeee-na-myetallurgeechyeskeeh-proeezvodstvah/?ysclid-Inshtjh0z2208396380 (дата обращения: 22.11.2023). 

3ГОСТ 12.0.03-2015. Система стандартов безопасности труда. Опасные и вредные производственные факторы. Электронный фонд 
правовых и нормативно-технических документов. URL: https://docs.cntd.ru/document/1200136071 (дата обращения: 21.11.2023). 
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На графиках рис. 1 достаточно выражена синхронность пиков, а также минимальная аварийность B 
последние часы работы. Значит, время работы оператора — существенный фактор возникновения несчастного 
случая на производстве. 

В [3] экспериментально установлена необходимость учитывать хронотип оператора при определении 
профессионально важных качеств. 

Для предприятий металлургической отрасли методы математического моделирования эффективны при 
анализе ситуаций, которые: 

— связаны с промышленной безопасностью; 

— могут возникнуть в сложной технической системе «человек — машина»; 

— формируются и развиваются в различных условиях внутренней и внешней среды. 

Математические модели, описывающие психофизиологическое состояние оператора, разрабатывались в 
различных математических постановках. При этом задействовали инструментарий измерения первичных 
показателей. В.Г. Абашин выяснил, как использование биометрических технологий связано с 
психофизиологическим состоянием оператора автоматического рабочего места, его работоспособностью, 
снижением количества брака и аварий технологического процесса по вине человека [4]. Работоспособность 
оператора моделируется по клавиатурному почерку. С этой целью применяют нечеткие множества как основу 
адаптивных моделей — биометрической и мультибиометрической [5]. 

Авторы [6] используют в математической модели биометрические признаки (голос, клавиатурный почерк и 
характер работы с компьютерной мышью). Это позволяет судить о психофизиологическом состоянии 
оператора: норма, усталость, опьянение, возбуждение, расслабление (засыпание). Модели базируются на 
стратегии Байеса, а также нейросетевом подходе и позволяют оценить уровень ПФС оператора и 
спрогнозировать его способность к выполнению текущих операций. 

В [7] разработаны имитационные математические модели взаимосвязи факторов психофизиологического 
состояния оператора. Авторы задействовали рекуррентные соотношения в дискретной модели и системы 
дифференциальных уравнений в непрерывной модели, что позволяет ввести в рассмотрение фактор времени. В 
этом случае параметры модели определяются по: 

— результатам тестирования конкретного работника; 

— его способности к определенному виду деятельности; 

— итогам экспериментальных исследований рабочего места и функциональных действий оператора. 

Эта часть моделирования меньше других формализуется и в наибольшей степени влияет на возможности 
использования модели для корректировки реальных действий оператора. В настоящее время проблема не 
решена. Таким образом, сохраняют актуальность вопросы определения параметров математических моделей, 
которые могут корректно описать взаимосвязи различных факторов психофизиологического состояния 
оператора. Цель данной статьи — разработка методики определения параметров таких моделей. Научные 
изыскания в данной сфере будут способствовать выявлению показателей состояния оператора 
металлургического оборудования, которые влияют на его функциональные характеристики и могут стать 
причиной серьезных ошибок с негативными последствиями. 

Материалы и методы. Для оценки психофизиологического состояния оператора выбрали следующие 
показатели: 

— работоспособность (возможность целенаправленно выполнять работу определенного количества и 
качества за обозначенное время); 

— утомляемость (снижение работоспособности с нарушением координации движений, падением 
концентрации внимания и точности решений [8]); 

— ошибаемость (оценивается количеством ошибочных действий). 

Утомление определяется по субъективным и объективным признакам [9]. Субъективный уровень утомления 
устанавливали по шкале оценки усталости“ (fatigue assessment scale, РАЗ) [10]. Ее разработала группа 
нидерландских ученых под руководством X. Михильсен (H.J. Michielsen). Это 10 вопросов о ежедневном 
самочувствии респондента. 

В качестве источника объективных параметров в данной работе использовали тест цифровой корректурной 
пробы (ЦКП) По его результатам вычислялись экспериментальные значения показателей 
психофизиологического состояния оператора. Числовые значения параметров модели определяются по 
названным показателям для конкретного респондента в совокупности с используемой моделью. 


* Шкала оценки усталости: Fatigue Assessment Scale (FAS). URL: https://wasog.org/dynamic/media/78/documents/Questionairres/fas rus anon.html (дата 
обращения: 21.11.2023). 
5 Цифровая корректурная проба. URL: https://metodorf.ru/tests/korrekt/korrektchis.php (дата обращения: 18.08.2023). 
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Задача подбора параметров математической модели взаимосвязи различных психофизиологических 
факторов оператора решается в два этапа. 

1-й — экспериментальный. Участников тестируют, собирают данные об условиях работы оператора и его 
личностных качествах, которые могут влиять на исследуемые параметры трудовой деятельности. 

2-й — расчетный. Используется база данных, собранных на 1-м этапе. С учетом выбранной математической 
модели оцениваются ее параметры для каждого конкретного работника. При достаточном объеме 
статистических данных в широком диапазоне варьируемых параметров строятся усредненные модели, которые 
можно применить для прогнозирования изменения психофизиологических показателей работы оператора в 
различных условиях и с учетом временного фактора. 

В данном исследовании для электронного тестирования использовали программу, представленную 
интернет-ресурсом «Интерактивный портал — книга методик саморазвития и достижения успеха». 

При выполнении ЦКП испытуемый строка за строкой просматривал сгенерированный программой числовой 
массив и зачеркивал цифры, указанные в задании. На выполнение теста давали 3 минуты. Затем программа 
выдала описанные ниже результаты. 

1. Основные первичные показатели: 

— время выполнения теста f; 

— общее количество просмотренных цифр до последней выбранной цифры № 

— общее количество просмотренных строк С; 

— общее количество цифр, которые необходимо вычеркнуть и; 

— общее количество вычеркнутых цифр М; 

— число верных ответов 5; 

— число пропущенных цифр Р; 

— число ошибочно выбранных цифр О. 

2. Расчетные показатели, характеризующие: 

— скорость (производительность) внимания А; 

— точность работы 7 (в трех вариантах); 

— умственную продуктивность Е; 

— умственную работоспособность Ан; 

— концентрацию внимания К; 

— устойчивость концентрации внимания Ky; 

— объем зрительной информации V; 

— скорость переработки информации О. 

Для решения поставленной задачи отобрали показатели, на основе которых оценивались 
работоспособность X, утомляемость Уи ошибаемость Z оператора по формулам (1—4). 


1. Умственная работоспособность А, (по цифровой корректурной пробе”, единица измерения — знаков в секунду): 


A, -(N /t)x((M -(O& P))/n). а) 


Для приведения показателя к нормализованной шкале [0, 1] в безразмерных единицах авторы предлагают 
формулу (2). Результатом будет i-roe значение показателя работоспособности X (соответствует дискретному 
условному времени ?): 


х= Аш {тах { А, } р (2) 


где max {Áu} — максимально возможное значение умственной работоспособности по данному виду теста 
независимо от тестируемого работника. 

Рассмотрим стандартный тест ЦКП с 1600 цифрами. Допустим, респондент за 180 секунд полностью и 
правильно выполнит тест. В этом случае величина тах{А,} составит 8,889 знаков в секунду. Если же 
тестируемый выполнит тест менее чем за 180 с, значение X при расчете по формуле (2) может оказаться 
больше 1. Тогда для Х принимается значение 1. 

2. Утомляемость У оценивается по показателю К (концентрация внимания): 


Y -(1-K /100), (3) 


где К = (М-0) п 100 — коэффициент, характеризующий концентрацию внимания, %. 


* Интерактивный портал — книга методик саморазвития и достижения успеха. URL: https://metodorf.ru/ (дата обращения: 18.10.2023). 
7 Цифровая корректурная проба. URL: https://metodorf.ru/tests/korrekt/korrektchis.php (дата обращения: 18.10.2023). 
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Данный показатель У принимает безразмерные значения в интервале [0, 1]. 
3. Ошибаемость Z предлагается оценить по формуле: 


Z=(0+P)/(M +P). (4) 

Данный показатель также принимает значения B интервале [0, 1]. 
В дальнейшем будем характеризовать Х, У и 7 как унифицированные количественные индикаторы 
психофизиологического состояния оператора, связанные с его основной производственной деятельностью. Для 


трактовки показателей Х, Уи 7 разработана шкала, представленная в таблице 1. 


Таблица 1 
Интерпретация значений показателей в качественной форме 
Показатель Интервалы значений показателей / качественная характеристика 
[0-0,20] [0,20—0,40] [0,40—0,60] [0,60—0,80] [0,80—1,00] 
Работоспособность X 
низкая ниже средней средняя хорошая высокая 
[0—0,20] [0,20—0,40] [0,40—0,60] [0,60—0,80] [0,80—1,00] 
Утомляемость У 
низкая ниже средней средняя выше средней высокая 
[0-0,01] [0,01—0,05] [0,05—0,10] [0,10-0,20] [0,20-1,00] 
Ошибаемость 7 
незначительная заметная существенная значительная высокая 


Показатели X, Уи Z — безразмерные, унифицированные со значениями в интервале [0, 1]. 

На втором этапе для подбора параметров математической модели есть два способа. 

1. При помощи модуля «Поиск решения» программы Excel создается сценарий решения системы для 
некоторых начальных параметров. В его основе — система рекуррентных уравнений математической модели в 
дискретной форме. Затем «Поиск решения» выдает оптимальные значения параметров модели для целевой 
функции, равной сумме квадратов отклонений расчетных значений показателей X, Уи Z от экспериментальных. 

2.В программе Anylogic*, которая позволяет в автоматическом режиме подбирать параметры 
рассматриваемых математических моделей в дискретной форме рекуррентных соотношений и в виде системы 
дифференциальных уравнений. —Anylogic целесообразно использовать после анализа результатов 
предварительной оценки параметров модели с помощью модуля «Поиск решения» Excel. 

Для апробации предлагаемой методики отобрали 10 человек — студентов и сотрудников Донбасского 
государственного технического университета. Они работали как операторы на учебных тренажерах с 
автоматическим и полуавтоматическим управлением. В начале и в конце рабочего дня оценивался 
субъективный уровень утомления по шкале ЕА$. В начале каждого часа, с 9.00 до 16.00, выполнялись тесты 
ЦКП. Данные обрабатывали по формулам (1-4) для каждого респондента. Результаты тестирования и их 
унификацию (X, Уи Z) использовали на втором этапе для определения параметров математической модели. 


В данной работе применили модель в рекуррентной форме [7], которая описывается системой уравнений: 


X. х. 

Хи=х, +a —— (1-7 2) - byx;y; = hx;z;, 

+1 dus ^e ( x! Ху 
Yi Yi 

"LSU : (1-59 by; 5 

Xia = У ? dy е ( m. 2i y. (5) 
2. 7. 

Z- — Z; +a; ————-:(1- —)-bhiy;z.. 

1+1 i 3 32, + © ( m ЗУ; i 


Решением системы (5) являются три сопряженных временных ряда длиной т индикаторов 
X = (x, 1=0, 1, 2, ..., т}, Y= {у, 1=0, 1, 2,..,т),  Z-71z,i-0,1,2,..., т). Индекс: — переменная, 
характеризующая дискретное время в системе. Параметры aj bj К, h, 4, с; ( = 1,2, 3) определяются в 


результате решения оптимизационной задачи по исходным данным тестирования. 


$ AnyLogic: имитационное моделирование для бизнеса. URL: https://www.anylogic.ru/ (дата обращения: 12.11.2023). 
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Математическая модель оптимизационной задачи описывается системой ограничений (6) и целевой 
функцией (7). 
8,(aj b, kj dj cj) e GP х)? SA, 
8s(aj bk hid; e; ) e (у-у? < 4, 
б, (аЬ. су) (20-20? < А, 
a by, kj, hd ;,c; 2 0, o 
хР 20, y? 2 0,z? 20, 
хр <1,ур <1,z? <1, 


i=1,...,m, 
j=1,2,3. 


V(aj bj kj, h dye; ) - У, +У 5, +} б > тм. (7) 
і=1 m izl 


Здесь индексом ф помечены фактические значения индикаторов X, Уи Z, найденные в результате обработки 
данных тестирования оператора, р — расчетные значения Х, Уи 7, определяемые из решения системы (5) и 
зависящие OT искомых в данной задаче оптимизации переменных ау, Бу, К, h, dj, с; (ў = 1, 2, 3). Ha варьируемые 
переменные aj, bj, Ку, А, dj, с; ограничения задаются по максимально допустимому значению À квадратов 
отклонений (Öxi, бу, 04) расчетных значений X, Уи Z от фактических для всех т значений рассматриваемых 
временных рядов. Целевая функция V(a;, Ру, kj, h, dj, cj) равна сумме квадратов отклонений (Öxi, ду, Өш). 

Для решения оптимизационной задачи применялся метод обобщенного приведенного градиента (ОПГ). 

Результаты исследования. 1-й этап — экспериментальный. Результаты тестирования участников 
эксперимента (респонденты №1-№10), преобразованные к унифицированным показателям X, Уи Z, представлены 
на диаграммах размаха (рис. 2-4). В виде точки показан центр распределения (медиана), прямоугольник 
обозначает границы варьирования (квартили 25 % — 75 %), «усики» — нижнюю и верхнюю границу значений 


показателя (min — max). Диаграммы выполнены с помощью программного комплекса «Statistica». 
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Рис. 2. Показатели работоспособности операторов 
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Рис. 3. Показатели утомляемости операторов 
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Рис. 4. Показатели ошибаемости операторов 


Анализ рис. 2-4 позволил выделить респондентов с лучшими и худшими оценками их деятельности 
(таблица 2). Учитывались уровни показателей Х, Уи 7, соотнесенные с градациями таблицы 1, а также квартили 


и интервалы варьирования показателей. 
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Таблица 2 
Ранжирование респондентов по результатам тестирования 
E Варианты (Ne респондента) 
Критерий оценки 
наихудшие наилучшие 

Работоспособность Х (чем больше, тем лучше) 8,1,4,6 7,9,3 
Утомляемость Y (чем больше, тем хуже) 8, 1, 4,6 7, 10 
Ошибаемость Z (чем больше, тем хуже) 4,8, 1 9, 10,7 
Интегрально (с учетом Х, У, 7) 8,4, 1 7, 10, 9 


Графики на рис. 5 отражают динамику работоспособности для респондентов №7 (пунктирная линия, 
лучший результат) и № 8 (сплошная линия, худший результат). 
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4 07 e---9—79. --- > 
o 9 `~ > 
xx 2 r 
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9 T 27 > 
5 А a 
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Рис. 5. Динамика работоспособности операторов 


Для объяснения достаточно большой разницы данных респондентов 7 и 8 учитывали их возраст, пол, род 
занятий, самочувствие и др. Наиболее существенным оказался фактор хронотипа. Классификация, принятая в 
1970 году, подразумевает три категории людей с отличительными поведенческими признаками и генетической 
разницей биоритмов. Это так называемые «жаворонки», «совы» и «голуби» [11]. Респондент № 8по 
хронотипу — «жаворонок», №7 — «сова». В результате работоспособность № 8 снижается к концу дня, 
а №7 — повышается (рис. 5). Хронофизиологические особенности операторов проявляются также в динамике 


их утомляемости и ошибаемости. 


2-й этап — расчетный. Получены параметры модели (5), рассчитанные для каждого участника 
эксперимента. С этой целью решили оптимизационную задачу (6)-(7). Диапазоны изменения параметров 


приведены в таблице 3. 


Таблица 3 


Диапазоны параметров математической модели (5) для группы участников эксперимента 


Показатель Диапазоны параметров Показатель Диапазоны параметров 

а, 0,0004—0,0300 bs 0,0001—10,1400 
bi 0—0,9700 Кз 0,0020—0,0700 
kı 0,5900—1,0000 Сі 0,0010-0,2000 
h 0—1,8100 с? 0,0001—0,0010 
az 0—0,0010 Сз 0,0100—0,3300 
b» 0,0070—0,2600 di = 0 

k2 0,4200—1,0000 d» = 0 

аз 0,0005—0,0600 d3 0,0400—0,9800 


В таблице 4 представлены результаты оценки качества математических моделей динамики изменения 
показателей X, Уи Z, оцениваемых по среднеквадратической ошибке модели и средней относительной ошибке 


для всех испытуемых. 
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Таблица 4 
Оценка качества математических моделей 


№ респондента 1 2 3 4 5 6 7 8 9 10 


Средняя относительная ошибка =, % | 14,27 | 4,54 | 20,93 | 12,80 | 7,28 | 12,00 | 9,66 | 15,03 | 12,23 | 9,70 


Среднеквадратическая ошибка, ед. | 0,182 | 0,107 | 0,253 | 0,267 | 0,158 | 0,247 | 0,271 | 0,256 | 0,192 | 0,166 


Для разных респондентов показатели качества модели варьируются в широком диапазоне. Отличное 
качество (= < 5%) получено в одном случае, хорошее (= < 10 %) — в трех случаях, удовлетворительное 
(є < 15%) — в четырех. Для двух случаев (№3 и №8) не удалось удовлетворительно решить задачу 
оптимального подбора параметров математической модели. Это говорит о том, что не все влияющие факторы 
учтены либо модель (5) в некоторых случаях не работает. 

Согласно результатам анализа алгоритма решения поставленной задачи, ряд упрощений улучшит 
сходимость результатов. Можно, например: 

— уменьшить на 3 ед. число параметров системы (5), сократив числители и знаменатели вторых слагаемых 
правых частей системы на C1, C2, сз соответственно; 

— удалить из системы (6) ее первое, второе и третье ограничения, переведя их в статус наблюдаемых ограничений. 

В итоге алгоритм метода ОПГ будет работать лучше. 

На рис. 6 и 7 приведены сравнительные диаграммы оценки качества построенных математических моделей 
в двух вариантах: 

— ] — исходная модель (5); 

— 2 — модифицированная модель. 

Как видно из диаграмм, решение задачи подбора параметров математической модели динамики показателей 
ПФС оператора значительно улучшено при использовании модели 2 для всех случаев, кроме № 8. 
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Рис. 6. Сравнение математических моделей по средней квадратической ошибке 
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Рис. 7. Сравнение математических моделей по средней относительной ошибке 
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Обсуждение и заключение. Предложенный подход позволяет получить и прогнозировать для каждого 
индивидуума динамические профили его психофизиологических характеристик, составить математические 
модели взаимосвязей. Целесообразно увеличить объем выборки, уменьшить шаг дискретного времени и 
выполнить исследования для различных условий работы, включающие технологические, климатические, 
экологические, психологические и прочие факторы, которые могут влиять на эффективность работника. Это 
следует сделать для расширения функциональных возможностей моделей, их применения в реальных условиях 
работы оператора металлургического оборудования. 

Обнаруженный эффект влияния хронотипа человека на его ПФС требует более основательных 
исследований. В статье [12], посвященной поиску компонентов циркадных часов у людей, на больших 
статистических выборках показано, что хронотип зависит от многих факторов (пол, возраст, график работы и 
др.). Все это необходимо диагностировать и учитывать в математических моделях оценки ПФС операторов. 

Интересно также дополнить модель компонентами культуры производственной безопасности [13], которые 
можно привести к количественному виду и использовать как поправочные коэффициенты. 

В дальнейшем предполагается усовершенствовать модель дополнительными переменными. Для этого нужно 
зафиксировать психофизиологические показатели и определять местонахождения оператора в режиме 
реального времени. 

Представленную в работе методику можно взять за основу при решении описанных ниже задач. 

— Составление и анализ динамического профиля работника при приеме на работу в качестве оператора 
машины, агрегата или устройства, где важны такие характеристики ПФС, как работоспособность, утомляемость 
и ошибаемость. Эта задача решается для конкретного человека в определенных условиях производства, что 
позволяет подобрать оптимальный режим работы и отдыха оператора, способствует сохранению здоровья 
специалиста и повышению уровня промышленной безопасности предприятия. 

— Формирование и ведение специализированных баз статистических данных, включающих характеристики 
рабочих мест и сотрудников, в том числе их психофизиологическое состояние. На основе собранного 
статистического материала и систем мониторинга ПФС оператора в режиме реального времени можно строить 
прогнозные модели для предупреждения нештатных и аварийных ситуаций на металлургических предприятиях. 
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